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Parameterized temporal feature analysis 



The invention relates to classifying at least one audio signal into at least one 

audio class. 

5 Developments in Internet and broadcast technology enable users to enjoy large 

amounts of multimedia content. With this rapidly increasing amount of data, users require 
automatic methods to filter, process and store incoming data. Some of these functions will be 
aided by attached metadata that provides information about the content. However, due to the 
fact that metadata is sometimes not provided, not precise enough, or even inaccurate, and 

10 because local processing power has increased tremendously, interest in local automatic 
multimedia analysis has increased. The multimedia analysis may comprise an automatic 
classification of an audio signal. In the automatic classification of an audio signal, low-level 
signal features are mapped to a semantic meaning, i.e. a classification of the analyzed audio 
content. By way of example and not limitation, the classification may be a discrimination 

15 between music, speech, background noise, or silence. Also other classifications are possible, 
such as musical genre classification, automatic detection of particular audio highlights or 
automatic speaker recognition. The classification of audio usually comprises two stages. The 
first stage analyzes the incoming waveform and extracts at least one audio feature that gives 
information about a predetermined property of the audio signal. The feature extraction 

20 process usually involves a large information reduction. The second stage performs a 
classification based on the extracted audio features. 

E. Wold et al. presented a number of audio features that can be used for an 
audio classification in "Content-based classification, search and retrieval of audio" in IEEE 
Multimedia, Fall: 27-36, 1996. These audio features are loudness, pitch, brightness and 

25 bandwidth of an audio signal. The audio features may vary over time, which results in so 
called trajectories of the audio features. In order to obtain information about the feature 
trajectories, i.e. the temporal fluctuation of an audio feature, a number of further audio 
features are introduced. These further audio features comprise the average value of an audio 
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feature over a feature trajectory, the variance of an audio feature over a feature trajectory, and 
the autocorrelation of an audio feature over a feature trajectory at a small lag. 



5 It is an object of the present invention to obtain information about the temporal 

fluctuation of an audio feature in an advantageous manner. To this end, the invention 
provides a method, a system, a music system, a multi-media system and a medium as 
specified in the independent claims. Advantageous embodiments are defined in the dependent 
claims. 

10 According to an aspect of the invention, to obtain information about the 

temporal behavior of an audio feature, a further audio feature is introduced which is based on 
a frequency analysis performed on a set of values of said audio feature at different time 
instances. The invention is based on the insight that, although the average and the variance of 
an audio feature over time does provide information about the temporal fluctuation of an 

15 audio feature, it does not provide any information about the temporal fluctuation velocity of 
the audio feature. Moreover, the average and the variance of an audio feature are usually 
correlated. For example, if a signal is scaled by a factor of two, both the average and the 
standard deviation of the short-term energy scale with the same factor. Most classification 
algorithms work more efficiently if the analyzed audio features are uncorrelated. Next, the 

20 autocorrelation of an audio feature introduced in Wold et al. may be a measure of whether or 
not an audio feature is changing over time and whether or not it is periodic. However, the 
autocorrelation does not give a detailed description of the temporal behavior of an audio 
feature. The autocorrelation may give an indication of how fast an audio feature is varying in 
time, but this indication is averaged over the whole signal. Therefore, using the 

25 autocorrelation of an audio feature is only giving limited information about the temporal 
fluctuation of an audio feature. Using a further audio feature according to the invention 
solves at least one of the disadvantages stated above. 

According to a further aspect of the invention, the audio feature that serves as 
an input for the frequency analysis may be at least one audio feature that is known in the art. 

30 By way of example and not limitation, an audio feature can be chosen from a plurality of 
audio features, such as a root-mean-square (RMS) level, a spectral centroid, a bandwidth, a 
zero-crossing rate, a spectral roll-off frequency, a band energy ratio, a delta spectrum 
magnitude, a pitch and a pitch strength. These audio features are commonly used features that 
are known in the art. An advantage of using these audio features is that it is relatively simple 
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to calculate them which is advantageous for the required computational load. A further 
possibility to choose an audio feature is to use at least one mel-frequency cepstral coefficient 
(MFCC). MFCC coefficients represent a parameterized description of the amplitude 
spectrum of an audio signal. An MFCC coefficient is used in audio classification algorithms 
5 due to its compactness, i.e. MFCC coefficients are able to represent the spectral envelope 
with only a few parameters. Furthermore, the MFCC coefficients are approximately 
uncorrelated for speech signals and music. Also, except for the zeroth MFCC coefficient, 
which is a function of the overall signal level, the remaining coefficients do not depend on 
the input level, i.e. they are gain independent. A still further possibility to choose an audio 

10 feature, is to use common known psycho-acoustic features. By way of example and not 

limitation, these features can be the loudness and sharpness of an audio signal. Loudness is 
the sensation of intensity and sharpness is a perception related to the spectral density and the 
relative strength of high-frequency energy. Choosing these features for obtaining the further 
audio features may be advantageous as the psycho-acoustic features are related to a human's 

1 5 perception of audio. 

In an embodiment of the invention, in order to derive the further audio feature, 
an average (DC) value is calculated of a set of values of an audio feature at different time 
instances, at least one frequency band is defined, the amount of energy within said frequency 
band is calculated from said frequency analysis; and said further audio feature is defined as 

20 said amount of energy in dependence on said average (DC) value. An advantage of using a 
frequency band is that this frequency band may be chosen to correspond to a specific 
perceptual phenomena that may be important for audio classification. For example, speech 
signals contain prominent envelope modulations in the range of 3-15 Hz, which range 
corresponds to the syllabic rate. Other audio signals, such as music audio signals, have 

25 relatively fewer modulations in this range. Therefore, if speech audio signals need to be 

classified, it may be advantageous to use a further audio feature representing the amount of 
envelope modulation in the range of 3-15 Hz. Furthermore, envelope modulations in the 20- 
150 Hz range are perceived as roughness, i.e. musical dissonance. Therefore, in order to 
distinguish dissonant or rough sounds from consonant or smooth sounds, it may be 

30 advantageous to use a further audio feature representing the amount of envelope modulation 
in the range of 20-150 Hz. Next, envelope modulations at very low frequencies, for example 
in the range of 1-2 Hz are perceived as changes in loudness. Therefore, in order to distinguish 
sounds with different rates of loudness changes, it is advantageous to use a further audio 
feature representing the amount of envelope modulation in the range of 1-2 Hz. Also, musical 
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tempo information is represented in the range of 1-2 Hz. It is noted that above mentioned 
frequency bands are given by way of example and not limitation. Other frequency bands may 
be chosen without departing from the scope of the invention. It is further noted that the 
frequency bands may be overlapping and may vary in time possibly in dependence on the 
5 audio signal, processing results, other external or internal parameters, or a combination 
thereof. 

In a further embodiment of the invention the further audio feature is 
determined by deriving at least one coefficient by performing a discrete cosine 
transformation (DCT) on the result of said frequency analysis. An advantage of using at least 

10 one DCT coefficient is that they are independent of the signal level. Furthermore, DCT 

coefficients may be highly uncorrelated which may be advantageous for audio classification. 
Also, with an increasing number of DCT coefficients, more details of the result of the 
frequency analysis are covered. In that manner, one can choose the detail level in 
combination with the resulting processing load. 

15 The aforementioned and other aspects of the invention will be apparent from 

and elucidated with reference to the embodiments described hereinafter. 



In the drawings: 

20 Fig. 1 shows a block diagram representing an embodiment of the method of 

classifying an audio signal according to the invention. 

Fig. 2 shows an embodiment of a music system according to the invention. 
Fig. 3 shows an embodiment of a multi-media system according to the 

invention. 

25 



Fig. 1 shows a block diagram 1 representing an embodiment of the process of 
classifying an audio signal according to the invention. As an input of the process, an audio 
signal A is used. The audio signal A may be a frame of PCM samples xfnj of length TV/. The 
30 audio signal A serves as an input for a feature extraction step 10. In the feature extraction step 
10, at least one predetermined audio feature F is determined. A possibility is that the 
extracted audio feature F is at least one from the following audio features, i.e. a root-mean- 
square (RMS) level, a spectral centroid, a bandwidth, a zero-crossing rate, a spectral roll-off 
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frequency, a band energy ratio, a delta spectrum magnitude, a pitch and a pitch strength. The 
RMS level of an audio frame of length N can be calculated as: 



A spectral centroid is based on a power spectrum P/fc/ of the audio signal A. 
5 The power spectrum P[k] may be obtained by an FFT operation: 

P W = |^2*»]exp(2^/J\o| (2) 

where k is the number of the power spectrum bin, which relates to the frequency /according 
to 

/l*3-*k (3) 

10 wherej£ is the sampling rate of the input signal. The spectral centroid S/may be defined as 
the center of mass of the power spectrum P[k] : 

The bandwidth Bj of the power spectrum P[k] may be defined by: 



'i — stm — (5> 

15 A zero crossing rate i? z may be defined as the number of zero crossings of an 

audio signal A occurring in a predetermined time frame, for example the number of zero 
crossings per second. The spectral roll-off frequency f r may be defined as the frequency for 
which the energy below that frequency is a predetermined proportion p (0 <p < 1) of the total 
signal energy: 



20 f r = /[maxf £p[*] < p^F[k]\ 

i ^Jt=0 *=o J 



(6) 



The band-energy ratio B r may be defined as the relative amount of energy 
present in a predetermined frequency range fi-f 2 Hz: 
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The delta spectrum magnitude ft is a correlate of the change in the spectrum. If 
two subsequent time-frames have (normalized) power spectra Pi[k] and P/+;/5fc/, then the delta 
spectrum magnitude may be defined by: 

/- = \^(4W}-4Kmf (8) 

5 The pitch 7* may be calculated by taking the maximum in the autocorrelation 

function within a limited range of delays. The pitch strength S may be defined as the height 
of the maximum peak in the normalized autocorrelation function corresponding to the pitch 
value. 

Next to extracting above mentioned audio features, the extracted audio feature 
10 F may also be at least one mel-frequency cepstral coefficient (MFCC). For determining a 

MFCC coefficient, for a given audio frame xfnj 9 with 0 < n < N-l, the power spectrum can be 
computed for example by taking a Fourier transform of xfnj 9 resulting mX[k]\ 

X W = mh[n]e^^ /N) (9) 

where hfnj denotes a temporal window. An example of such a window is a Hanning window 
15 which is known in the art. The amplitude spectrum \XfkJ\ ofX[k] is multiplied by a set of 
filter kernels. The center frequencies of these filters have a constant separation on a mel- 
frequency scale f m in dependence on the frequency / which may be defined by: 

/„=25951og 10 (l + //700) (10) 
The input spectrum is converted to a mel-frequency spectrum using a 
20 filterbank with k n triangularly-shaped filters G[k t with a spacing and a bandwidth that is^ 
linearly spaced on the mel-frequency scale. The mel-frequency cepstrum is then given by the 
logarithm of the inner product of the filter kernel and the amplitude spectrum: 

C[k„ ] = log 10 &\X[k]\G[k 9 kj) (11) 

In order to obtain the mel-frequency cepstrum coefficients (MFCC) c[n], the 
25 discrete cosine transform of the mel-frequency spectrum is computed: 

c[n] = Xq*Jcos[«(* M ~1/2)7C/K] (12) 
*« 

A further possibility to choose an extracted audio feature F is to use at least 
one psycho-acoustic (PA) audio feature, such as loudness or sharpness of an audio signal. An 
example of defining loudness is presented by Eberhard Zwicker et al. in "Psychoacoustics: 
30 Facts and Models", volume 22 of Springer series on information sciences, Springer- Verlag, 
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Berlin, 2 nd edition, 1999. An example of defining sharpness is given in "Sharpness as an 
attribute of the timbre of steady sounds" inAcustica, 30:159-172, 1974. A plurality of 
methods are known in the art are known to extract psycho-acoustic features, that may be 
chosen for obtaining the further audio feature according to the invention. 
5 The further audio feature is obtained in a further audio feature extraction step 

12. In the further feature extraction step 12, a set of values of an audio feature F at different 
time instances is defined. Possibly, sets of values may have a fixed size or may be variable. 
Also, sets of values may be overlapping. Next, a frequency analysis is performed on the set 
of values F(t) of an audio feature F in time t=[0..T], which results in a feature log power 
10 spectrum \F(f)\ 2 : 

2 



fiff = 



jF(t)exp(-27qft)dt 



(13) 



In order to parameterize the feature log power spectrum, the frequency axis of 
the power spectrum is summarized by summing the energy within at least one predetermined 
frequency band. The summing may also be a weighted sum of energies, for example in 

15 dependence on frequency or in dependence on the energy itself. By way of example and not 
limitation, the predetermined frequency band can be at least one of the frequency bands 1-2 
Hz, 3-15 Hz, and 20-150 Hz. The 1-2 Hz frequency band may be preferable to distinguish 
sounds with different rates of loudness changes as envelope modulations at very low 
frequencies are perceived as changes in loudness. Also musical tempo information is 

20 available from this frequency range. The 3-15 Hz frequency band may be preferable for 

classifying speech signals which contain prominent envelope modulations in the range of 3- 
15 Hz, which range corresponds to the syllabic rate. Other audio signals, such as music audio 
signals, have relatively fewer modulations in this range. The 20-150 Hz frequency band may 
be preferable to distinguish dissonant or rough sounds from consonant or smooth sounds as 

25 envelope modulations in the 20-150 Hz range are perceived as roughness, i.e. musical 

dissonance. Finally, the amount of energy within a predetermined frequency band may be 
divided by the average (DC) of subsequent values of the audio feature to yield a relative 
modulation depth. The average may be obtained by evaluating the 0 Hz energy in the feature 
power spectrum \F\. The result of this calculation is a further audio feature F mo d that can be 

30 used for classifying an audio signal. Another method to parameterize the feature log power 
spectrum is to transform the log-power spectrum \F(f) \ 2 into at least one coefficient C(m) 
using a discrete cosine transformation (DCT): 
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cto-UiL^\Jpt£* (14) 

i { frfa ) \F(of 
in which f a and^ denotes the start and end frequency of a frequency band of interest. 
Usually, the upper frequency f b is half the sampling frequency of£. Now, the coefficient 
C(m) is used as a further audio feature F mod . C(0) denotes the total amount of modulations 
5 averages on a log scale, hence C(0) is related to the overall modulation depth. Due to the 
division of \F(x)\ 2 by \F(0)\ 2 the modulation depth parameter is independent of the signal 
level. Furthermore, it is noted that DCT coefficients may be highly uncorrelated which may 
be advantageous for audio classification. Also, it is noted that with an increasing number m 
of coefficients C(m), more details of the feature log-power spectrum \F(f)\ 2 are covered. 

10 The classification of an audio signal is performed in the audio classification 

step 14. Various methods of classifying audio based on audio features are known in the art. 
Any of these methods may be used in the audio classification step 14, now using a further 
audio feature F mo & possibly in combination with an audio feature F derived in the feature 
extraction step 10. By way of example and not limitation, the audio classification step 14 

15 comprises at least one of multivariate Gaussian models, Gaussion mixture models, self- 
organizing maps, neural networks, k-nearest neighbor schemes and hidden Markov models. 
These audio classification methods are known in the art and it goes beyond the scope of the 
invention to discuss them in detail. The result of the classification step 14 is at least one audio 
class AC. It is noted that the classification of audio may also comprise similarity 

20 measurements, e.g. finding similar music. The classification may further comprise a 
segmentation procedure, e.g. identifying scene changes using audio feautures. 

Fig. 2 shows an embodiment of a music system 2 according to the invention. 
The music system 2 comprises a medium 22 comprising audio data. The medium 22 may be 
a removable medium such as a Compact Disk or may be incorporated in the music system 2, 

25 for example a hard-disk incorporated in the system. The music system may denote an audio 
receiver. In that case, the medium 22 is at another location and the audio data is broadcasted 
to the audio receiver, for example by radio, internet, or satellite. The medium 22 may also be 
an audio data generating device such as a microphone or a musical instrument. The audio 
data is transferred to an audio data processing device 24. The audio data processing device 24 

30 may comprise an audio amplifier and may further comprise audio processing means for 

adapting the audio data and control means for controlling an audio classification process. The 
resulting signal from the audio data processing device 24 is transferred to a sound generating 
device 26 such as a speaker. The music system 2 further comprises an audio classification 
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system 20. The audio classification system extracts features from the audio data according to 
the invention and classifies the audio data into at least one predetermined audio class. The 
audio data may be obtained from the audio data processing device 24 or may be obtained 
directly from the medium 22. The resulting audio class or classes are transferred to the audio 
5 data processing device 24 for further processing. For example, based on the audio class the 
audio data can be played or not, or based on the audio class information may be given to the 
user about the content it is listening to. The features used for classification may be available 
from a signal accompanying the audio data, the signal comprising at least one further audio 
feature which is obtained according to the invention. When the music system 2 denotes an 

10 audio receiver, the audio receiver may further comprise means to receive a signal comprising 
at least one further audio feature which is obtained according to the invention, which further 
audio feature may be used for a classification. 

Fig. 3 shows an embodiment of a multi-media system 3 according to the 
invention. The multi-media system 3 comprises a medium 32 comprising audio and video 

15 data. There may be different mediums containing either audio or video data. The medium 

may be a removable medium, or may be incorporated in the system 3 or may be available via 
broadcasting. The audio and video data is transferred to an audio and video processing device 
34. The audio and video processing device 34 comprises means to process the audio and 
video for playing the audio and video data on an audio and video playing device 36. The 

20 audio and video playing device 36 may for example be a television, a monitor or a projector. 
The audio and video playing device 36 may be incorporated in the system 3 but this is not 
necessary. The audio and video processing device 34 may further comprise control means to 
control a video analysis system 38. The video analysis system 38 is capable of analyzing the 
video data obtained either directly from the medium 32 or from the audio and video 

25 processing device 34. By way of example and not limitation, the video analysis system 38 
may comprise a video decoder for decoding video data that is available in compressed form, 
e.g. an MPEG decoder, a commercial block detector, a movie detector or a video 
classification system, which examples are known in the art. Results of the video analysis 
system 38 are transferred to the audio and video processing device 34. The audio and video 

30 processing device 34 may further comprise means to control an audio classification system 
20. The audio classification system 20 functions in the same way as in the music system 2. 
Furthermore, the audio and video processing device 34 may comprise means to combine the 
results from the video analysis system 38 and the audio classification system 20. For 
example, when a movie is detected by the video analysis system 38, the audio classification 
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system 20 may classify speech in the movie into a class of speech from a certain actor that 
plays in the movie. 

It should be noted that the above-mentioned embodiments illustrate rather than 
limit the invention, and that those skilled in the art will be able to design many alternative 
5 embodiments without departing from the scope of the appended claims. In the claims, any 
reference signs between parenthesis shall not be construed as limiting the claim. The word 
'compromising* does not exclude the presence of other elements or steps than those listed in 
a claim. The invention can be implemented by means of hardware comprising several distinct 
elements, and by means of a suitably programmed computer. In a device claim enumerating 
10 several means, several of these means can be embodied by one and the same item of 

hardware. The mere fact that certain measures are recited in mutually different dependent 
claims does not indicate that a combination of these measures cannot be used to advantage. 



